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THE SECULAR PERTURBATIONS OF TWO PLANETS MOVING IN 
THE SAME PLANE; WITH APPLICATION TO 
JUPITER AND SATURN. 


By Dr. G. W. Hiri, Washington, D. C. 


The solution of this problem, when we restrict ourselves to the first 
powers of the eccentricities, is as old as Lagrange, and is well known. Le- 
verrier, in going over this ground, attempted to include the effect of the terms 
of three dimensions with respect to eccentricities and inclinations. But when 
his method was applied to the four interior planets of the solar system it led 
to results that were nugatory. This method being that of successive approxi- 
mations, the expressions for the unknowns obtained in the simplest form of 
the investigation were substituted in the terms of three dimensions ; in conse- 
quence, he arrived at the same linear differential equations as before, but now 
augmented by known terms. His difficulty, in the case of the four interior 
planets, arose from the appearance in the results of integrating divisors which 
might receive very small, or even zero, values within the range of uncertainty 
of the values of the planetary masses. 

As far as the general question is concerned, no one has attempted to 
push the investigation further. Under these circumstances I have thought it 
might be well to treat as completely as we can the very simple case where we 
have only two planets executing their motions in the same plane. Although 
we see here at a glance that the problem is reducible to quadratures, yet this 
taken by itself does not constitute a practical solution. Some difficulties are 
encountered in deriving from the quadratures series suitable for calculating 
the values of the unknowns. These difficulties I have succeeded in surmount- 
ing by a process which would not suggest itself, I think, at first sight. 

In the application which I have made to the case of Jupiter and Saturn 
with neglected mutual inclination, I have carried the approximation to quan- 


* Annales de l’Observatoire de Paris, Tom. II, pp. 105-170 and pp. [38 51 
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tities of the fifth order, inclusive; and it is not difficult to see what must be 
done if it is desired to go further. 


The first thing to be done in this investigation is to find a proper develop- 
ment of the potential or perturbative function. Quantities belonging to the 
interior planet will be denoted by symbols without an accent, and those 
belonging to the exterior by symbols having an accent. Let, then, ~, 7, «@, 9; 
vw, and f denote severally the mass of the planet, the radius, the semi-axis 
major, the mean, eccentric, and true anomalies, while we denote the distance 
between the planets by —.. The potential function 2 is then given by the 
double definite integral 

0. 1 ree snl dg, 


4z* Jo Jo 


or, if the integration is accomplished with reference to the eccentric anomalies, 
by the double definite integral 


l nan PY mm : 
Q - , f du du’. 


ad 


These formule show that the potential function is proportional to the average 
value of the reciprocal of the distance when the mean anomalies are regarded 
as the independent variables, or to the average value of the product of the 
radii divided by the distance when the eccentric anomalies are the indepen- 
dent variables. As the eccentricities e and e and the lengitudes of the peri- 
helia ® and ©’ are the variable quantities whose forms as functions of the 
time we are seeking, it is plain they must be left indeterminate in the expres- 


sion we obtain for 2. Since can be expressed in terms of uw and w’ as a 
finite form, the second formula for 2 is to be preferred. 
If 7 be put for @ o’, the expression for in the case we treat, is 


ee 
r ° - r~ 13 
Fil —g - COS (7 7 7) = : 
7 . ‘ ide 


Thus, the expression for 2 becomes 
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If B denote the same function of : that Laplace's } is of a. the ratio ot the 
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mean distances, we may write 


1 + & 
7 ; 7* l : 
k 2 , COS (J J i) | 5) S Bcos7(f j i) 
y” a ” 2 ' e 
] ) + 
D 2 B,¢ ; ' l 
= denoting the base of natural logarithms. If we make <"*—1= s, and put 
1 ] l e* é 
7 . w 
, 2 l a " 
from the equations 
7 a(l € COS @%), 7 COs ay “7 (COs e), ? sin t a) l > sin ca 
it is easy to derive 
@ 
7? ay, (1 ws) 1 
esv¥—! : a 
l Ws 
Thus 
7 1 j y x BR S w / Ss w re , , 
, 7 j l Os l wos 


Seeking now an expression for /, in terms of s and s’, we have 
(1 24 COS ¢ a*)—3 Vi AD eib 4 


(we omit Laplace’s subscript 4, as it is unnecessary for the purposes of dis- 
tinction). We can regard 6” as an approximate value of /,, and the true 
value can be developed in a convergent series by Maclaurin’s Theorem, if the 
perihelion radius of the exterior planet always exceeds the aphelion radius of 
the interior; that is, if 
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The augmentation which @ receives is 


w 
7 (1 ws) | 1 
| WS 
7 7 . 7 
/ . w 
7) | l Ws) l . 
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‘Thus 
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' . ° ° . . ° . . rT : 
Substituting this value of 4; in the expression given above for — , and multi- 


plying the result by 


mm? mi w 
3 - Y) (1 Ws) 1 . 
add a é 
and employing the symbol to denote the operation of taking the coefficient 


of s's” in the development of a function of s and s’ in a series of integral 


powers and products of s and s’, we shall have 


, ~ ¢ » &a4 i—k i}, J 
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Let us put 


This quantity is then a function of e. Let £’;’' be the same function of e’ that 


i) is of e. Then we can write 
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This constitutes the infinite series to be employed in this investigation, 
and it remains only to study the properties of the functions of e denoted by 
FE}. By expanding the binomial factors involved in /;’ and performing the 


operation denoted by \, we shall get 


ry fe 1) (7 Pee /) 
EY 1) l ) 4 
ae | 


, tak Ey Oe —7-—  56—8 
wre & 1.2 G+HDG+%) | 


The series within the brackets is a case of the hypergeometric series 


~ 


1 “a. “(ad 1) 3(3 Rh «= (Gd lila 2) 3(3 1) (3 2 
r a* : : 
l.; 1.2.7(7 -1) 1.2.3.7(7—- DG - 2) 


treated by Gauss in a memoir entitled “ Disquisitiones generales circa seriem 


infinitam, etc.”* This series gives the value of /;’ in terms of 7 and w, but 


it may readily be transformed into another expressed in terms of e. Adopt- 


ing Gauss’s notation for this species of series 

} ’ . . . ° 
J 4 i ( / a, a Ll, w*). 
But from Gauss’s equation [100], p. 225 of the volume quoted, 
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In consequence 
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*See Gauss, Werke, Band III, p. 123 
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It is remarkable that when ¢ and / are integers the value of £;’ is equivalent 
to a rational function of the two quantities e and | 1 e*. For, when isa 
positive integer, the series first given terminates after a finite number of terms. 
The same thing occurs in the second series when 7 —/ is not negative. By 


Gauss’s equation [82], p. 209 of the volume quoted, 





‘ere ger es ae ; 9\—-2%=1 z,( 7-2 +2 7—2+1 . 
i; or 7 ” 7, 1. 2 fal e2) » frie > x 5 9 1.2 
From this it follows that 
(2 1) (2 — 2)... (2 —9) [« 9, —-%=-15,{7—2+2 7—24+1 . . 
. 1.2...j g ee Fe a ai 
a—1)...(2— 7) [« (1—e) — 1 (7—J7 L)(@ 4—2)le}? 
BR Re 2 1.(7+1) 2 
(2 ) Ll)... ) 4) ¢}* 
1.2.(7+1)(7+2) {2 ae |. 4 
which affords a finite expression for //’ when 7 is negative. It will be noticed . 
that /;’’ = 0, when 7, not zero, is not greater than /. ‘ 
In order that the symmetry of the expression for 2 may be seen, we will 
write the development of this quantity at length without the employment of 
the summatory signs ; 
0 ma | Lg ] 
$2 oa 1 . ‘ 
| db” ' 
7) Pi Adg A ad hat 
da q 
é l 7 hed fi ” 
Z - (2,08? —BZOOR'S EOr’® 
: )- dia - 3 
l 1°h° te a 0 
e-_ (E%EF’ ©  3F OF EOE’ — EO} 
2 ) da’ . 
) : 
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mm (Same expression as above except that 2, /, and) ' 
; cos 
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a  (£ now take 1 as the upper index instead of 0.) 

mm (Same expression, except that 4, 4, and £”) niall 
«@ «now take 2 as the upper index instead of 0.) aye 
mm (Same expression, except that 4, “, and ££") cos 37 


«@ (now take 3 as the upper index instead of 0.) 
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It may be noticed that the terms in £,° 4°"), L724", BSL, Bk 
BOERS, BSR, ete. can be omitted in writing the expression, as the latter 


a 
factors of these products vanish. However the symmetry is more apparent 
when they are retained. 

The following table exhibits the values of all the /"s required in develop- 
ing 2 to the terms of the sixth order, inclusive. They are expre ssed as fune- 
tions of e, and the finite form is given as perhaps more interesting than the 


development in ascending powers of e. 
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In the present 
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investigation it will be more convenient to make use of a 
in powers of I 1 | °iz=@, 


\ 2 


in terms of e the values 


By substituting in 


- f? 64. 


0) (1 — 02) 







OF TWO PLANETS MOVING IN THE SAME PLANI LS5 


making, for the sake of brevity, 7 7) kK, and earrving the deve lopme nt to 


terins of the sixth order, inclusive, we obtain 


he (ff L) (A 2)(h (A { ) 
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t Ach L)(/ 2) ; 
2 1.2.()7 L)( 2 
/ } »\ kik | } j 2 } ae! 
Sf Liu+) 2.4.6 ' 


Or, particularizing with respect to /, 
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If, in the expression for 2, we eall the function of the eccentricities 


;, =" 1)! LBS . Ss ;' ; 
which multiplies —— a j,i 0 the coefficient of cos jy, .W,’, and denot- 
Ud P 
} ing the characteristic of finite differences with respect to the variable 7, it will 


be seen that we have 
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Then the expressions for M, ean be derived by considering the preceding 
expressions, taken alternately with the positive and negative sign, as the sue- 
cessive differences of these functions with respect to the index 7; and it will 
be advantageous to apply the process separately to each power and product 


of 4and #7. The exhibition of this follows: 
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Coetticients of 4°. 
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We ean now write the explicit development of 2 as follows: 
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In order to have as few functions of « to deal with as possible, we vathe) 
together all the terms having the same powers of #4 and # as factors Also 
it will serve our purposes better to have the development of £7 in powers 
of cos } than in cosines of multiples of +7. For convenience in writing we 

d'h . F 
denote «' by (7,7). We then put 
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Then, neglecting the term which is independent of #, #, and 7 for the 


reason that it is useless for our purposes, we shall have 
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mn : - J 
A OFY A OPO Ag" 
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A, WH? cos y. 

In order to make an application of the method to the case of Jupi- 
ter and Saturn, we take from Runkle’s Tables of the Coefficients of the 
Perturbative Function the values of log (), 7) corresponding to the urgument 


log “ 9 736741 i. 


i. ) U). 4 l. 7) Z. ) o. 
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4 O.873T099 0.8712079 O.8S884960 0.9011986 
5 1.5571487 1.5610571 L.5658245 L.5ST98073 
6 2. 3402885 2.3412199 2.5462289 2.555596] 


Making use of these values, we obtain for this special case 
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The portion of the subject which treats of the integration of certain dif- 
ferential equations is now to be attended to. Denoting the mass of the sun 


by JM, and putting 
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the differential equations which determine the eccentricities and positions of 


the perihelia of the two planets are 


aG ALY do AL 
dt do ’ dt dG F 
dG dL da AL 
dt doa’ dt de 
But since 2 involves @ and «’ only through 7 = 0 — o’, we have 


dL AQ 
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Hence 


Gr (7 a constant 


is an integral of the problem. This integral equation may be more suitably 
expressed in terms of the variables # and # which we have before employed. 


Then A’ denoting an arbitrary constant, and denoting the constant quantities 


The value of A’ is ascertained by substituting in the left member of this 
equation for # and # the values they have at a definite epoch. We can 
how reduce the pumber of variables in the problem from four to three by 


adopting a variable » to replace f# and “#, such that 
“ 4, A sin dy, HW 1 A cos dy. 


Jy remains always in the first quadrant. Denoting the angles of the eccen- 


tricities by ¢ and ¢’, the eccentricities are determined by the formule 
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Making the substitutions in 2 necessary to make it involve » instead of # and 
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m 2 becomes, then, divisible by A’, and, in order to simplify, we 


Ai. Therefore, if we write « for cos v and put 


(] 
y\4 yg | 
(2 
(] - 
( Fi / 
») 
(4 
(] 
(“As A 
(4 
(] 
4 / | / 
(3 
4/4 l 
(7 
») ) 
_— { ~ 
4 AY ' 
( AA 
» » 
_ (- 
AK AY ‘ 
; ca 
_ ,4 
 « n 
AA 


n have 


R L, 





sin’ » cos’ 


SAME PLANI 


AKA 
Ci 2 
WHOA 
= 
hy? 
) 
ti 2 
LIA 


sin Y COS j 


/ 


/ 


A 


A 


A 


A 
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With this expression for /2 it is readily seen from the preceding differen- 


tial equations that the differential equation determining v is 


dy Ll dh 
dt sin vy dy ’ 
or 
da AR 
dt dy , 
Since 7? a constant is evidently an integral of the problem, we shall have 


Ake dy IR dy 


0). 
dy dat dy at 


W he nee is derived 
dy l ale 
dt sin vy dy — 


We still need an additional equation giving the value of some other function 


of w and Ww than Dy) ” . If we select wo wo we have 
dia @ ) AQ AL 7 
dt dG ddr” 


If A’ is kept evident in the expressions for the various /’s, so that the par- 


tial derivatives of them with respect to this quantity may be taken, we shall 


have 
dL d( KR) dk k Th dy Ld(ARY 1 dR 
adaG’ dkh ad G dy dG? » dh 9 y dy ° 
2 tan - 
di A(ANR\dkh k IR dy Ld(AR) 1 tans ® dk 
dG dk dy’ dy addi Z dh y. 2 dy : 
Whence 4 
J(a+-@) d (KR) cosy dk i 
dt dh sin » 4 
dl | A?) cos } dy q 
dh i Xs 


In making our numerical application we take the mean distance «’ as the 
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unit, when a@ becomes the same as «4 previously given, and assume for the 


masses the values 


l l 
a TL — , 
L047.879 3482.2 
These give 
log / L.5758667. log / 1.7708956. 


The values adopted for the eccentricities at the beginning of 1850 are 
0.04825801. 0.05606467 
These furnish the equations 


S.4778154 sin = 8.6728444| cos 5 


and the function /? becomes 


R 0.0005906465 0.00025439647 0. QOOOODOLIBTSO.2 
0.0000000193 94.22 
0.0003548741 0.00000629406.r 0.0000000873 Lr? | sin v cos 


0.000001 48778 0.00000004479.7! sin?» eos? yr 


/ 


0.000000006560 sin*» cos*;. 


The value of the constant in the integral equation 
a= C€ 


is ascertained by substituting in the expression for /? the values which » and ; 
have at a definite epoch, as 1850. C’ being determined, the equation /? — ¢ 
ean be solved, regarding sin » COS 7 aS the quantity whose value is to be 
obtained. This value can be supposed developed in powers of cos: vy 


and we write 


HT iD 


sin ¥ COS7 


0 


Dix D2" D,x* 


The readiest method of obtaining the //’s is by substituting the last expres- 
sion in # and then equating the resulting coefficients of each power of 2 to 
zero. We thus have a system of equations determining the //s. These can 
be solved by successive approximation. If (is allowed to appear as an inde- 
terminate in the expressions for the /’s, // can be partially differentiated with 
reference to this quantity. 
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We can now make // play the role of 72; for we have 


de dk dy ___ a ™ | il | 
; ; ane Yi are COS ‘ 4 
dt dy dt dar / ] l i 
Thus 
adil 
dt dy AC!’ 
dx Ac’ | ] 7) // 9? 
where the radical in the denominator must receive the sign of sin v3 for we 
have ‘ 
COS LV a sé 
sin» cos 7 // dD Dir Dix? Dr? 


/ 0 2 } ae 9 


Sin ? sin 7 


/ (1 — 2x — HH). 


If we suppose the orbits are always ellipses . cannot pass the limits & 
Thus . must oscillate between a maximum and a minimum value, while 
d/1/dC remains constantly of the same sign. The maximum and minimum 


values of . are evidently the two consecutive real roots of the equation in . 


1 xe* a * == %, 


which contain between them at any time the actual value of . Calling these 
roots a and b, we may write . 
l L* Hi * (a r) (xr b) (/, a 
a 
where ( is positive for all values of « lying between a and b; and when the a 
eccentricities are always small, the variation of ( is slight in comparison with of 
its magnitude. In the place of we can adopt a new variable, ¢', such that 
a b a b ; 
1 COS ¢ 
Z 2 
Then a 
ar dh MS 
} (a Ee. by) 
and the differential equation giving ¢’ in terms of 7 is Se 
dil 
dt dt’ 
] ae 
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To see how all this applies in the case of Jupiter and Saturn we assume 


the following values of the longitudes of the perihelia at the epoch 1850.0: 


. ” 11° 54’ 317.18, o 90° 6 5T”.55. 
The value of the constant ( being now determined, and the equation R ( 
modified in such a way that it becomes more suitable for solution, we have 
0.4021256 0.7168638x% 0.005545 Le* 0.0000546..% 
| 0.017 7360+ 0.0002460..7 | sin » cos ; 
. 0.0041924 0.0001262.c| sin?» cos*; 
} 


00000185 sin®y cos? y, 


When this equation is solved with reference to sin » COS 7 aS the unknown, we 


obtain 


// 0.4028046 0.7121389.~ O.0050141+" 0.QQ000050.2 


And when we ascertain what increment // receives from an infinitesimal inere- 
ment in the quantity C’, it results that 


Ail a rs es 
2827 .425 33.179. 0.005..7 
d¢ 
The equation l p> H ? 0. in this case, is 
O.8377485 0.5737057x 1.5031028.r 0.00714472" O.Q0000L8024 = 0. 


The consecutive real roots of this which contain between them the value of 2 
at 1850.0 are 





i a = 0.58038236, by O.9586738 
ee We derive from these the limiting values of », which are 
a 
Me 54° 31° 367.14 and 163° 28’ 14”.01, 
4 
i Thus, when 7 = 0°, the minimum e of Jupiter has place, which is 0.02752623 ; 
as also the maximum e’ of Saturn, which is 0.08362800. And, when ; 180 


the maximum e of Jupiter has place, and is 0.05944555; and the minimum « 
of Saturn, which is 0.01353514. 












































a 
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The remaining factor of the equation, two of whose roots we have just : 
obtained, is 
() = 1.5058180 — 0,0071522” — 0.0000180.°. 
Whence 
0.8149177 — 0.0019353.7 —- 0.00000202?. 
i ¢ 
Substituting, then, for 7 the expression 
r 0.1891751 — 0.7694987 cos ¢/, 4 
we get 
d = ~ yr ~ » : 
9304.1185 — 21.5662e — 0.05432 ’ 
de y 


Y308. 1802 16.5794 cos ¢ 0.0161 cos 2c. 
Integrating this, ¢ being the arbitrary constant, 
t-- ¢ —- 2808.1802¢ — 16.5794 sin ¢ — 0.0080 sin 2¢. 


Inverting this series and changing the numerical coefficients into seconds of 


are we get 


(f L9”.05825 (7 c) 1481°.57 sin |19°.05825 (¢ c) 


i 


6.04 sin 2 |19°.05825 (7 c)|\. 


From the value which ¢) must have at the epoch 1850.0, ¢ being counted 


thence, 


10” O5825e 977° ¥ 9’ .15. 
Also, we have 
COS » 0.189175] 0.7694987 cos , 
sin » COS 7 0.2679063 — 0.5494490 cos ¢ — 0.0029673 cos?’ 


0.0000023 cos*¢, 


sin» sin 7 0.9446898 0.0017265 cos ¢ 0.0000018 cos *¢/| sin h, 


These equations enable us to determine the eccentricities and difference of the 
longitudes of the perihelia at any given time. 
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It remains to find the longitudes of the perihelia themselves. We have 


da 1 , | ~a ee l - dy 
—3/( 4 - 
dt , 2 . dh 2» dt 
dt ...to4® 1-44 
~ ¢ ‘ 
dt . » in dh » dt- 
Or 
dia J (7) ] " dy l vil; 
dp 3 A dh ~? hy ; 
dia : (7) l - dy | rads 
oe / - / 
di 2 ’ dh 4 Lr 
Or 
pdH di Ma 
dja—tCh 1° dak br 1 
he 4 = l x fi? 
K’ dil ; AH Ile 
dia : Ct) | dh ha | / 
dr 2 } | | a //] 3 


Here A’ must be left indeterminate in the coetticients //,, 7), ete. of //, in order 


that we may get —- . In the next place, we derive 
dH dM Mh 
/ { w ( 7) l dah , dr ] / 
AS » i WV 
» dil dil Tha 
, A . (1 r) 
dla ( t) l dh da l 
dh 2 1 


When //, which is an infinite series in integral powers of .”, is divided by 1 

or 1 + 2, remainders independent of are left over which are equivalent to 
what //7 becomes when in it we make “ = 1 and + 1. These remainders 
we denote as //(1) and //(— 1). Then we may write 


H(1) T[X5x (t(—1)D,—t Da, + Drage Dy | 


dia ; C7) 1 ] r dh 
Ags > ] (/ 
H(-1\) ‘=P ,-dD,_ |. : 
/ SD ach — 9 \9P,-+-40Pi 0. -+- Pere — Par ore” 
d(a— Ct) 1 lia st ' dh : — ‘f 


dd 2 ] () 








arn eee © 
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The difference of these equations gives 


1M’) 1M—1) SP _ayp ap op |: 


_— +] 
J; 3] Bi 2 J r oc 


7h ] () 


Since + 
SiMmMce 7 


it follows that when the right member is expanded in an infinite series con- 


returns to the same value after ¢) has augmented by a circumference, 


taining, besides two terms in the form of fractions having 1 v and I ! 
as denominators, a set of terms proceeding according to cosines of multiples 
of ¢', the coefficient of the zero multiple of ¢’ must vanish. This is not imme- 
diately evident from the form of the expression. Hence I proceed to prove 


it to the degree of approximation we adopt. Let 


I 


E+ £e + Eo’ 
| (/ 5 
then, omitting the two terms in the form of fractions and having 1 r and 
| « for denominators, it will be perceived that we have 
(|: , ’ ’ ’ , ’ , 
“ DT: (D DA: DF, — (DE, D, F\ « 2D) bk: DE, \z*. 


Substituting for its value in terms of ¢/, if our proposition is true we ought 


to have 
DE, + (Dy + D,) E, + DB, — (DB, — Dk) * 
2D.E, DAK. E = ° ab | 0 
But if 
Q= M, + Mx + Mz? 
k= MA, 1 l y-iM,, } l VM, +? MM 
‘ . 6 S 


and M,, M, ; M,. a, and b are determined by the equations 


abd, = DI — 1, 
(a + b) M, — abM, 2DD,, 
M, —(a +b) M, 4+ abM,—1+ D? + 2D,D,, 
M, — (a+ b)H, = 2(DD, +- D,D;), 
3M. = D2 + 20,,D.,. 
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a By substituting the values of “#,, /,, and £, and multiplying by J/,!, our 
a equation becomes 
' p.* b DD k 3\* : b ib] V 
( “ ~ ) 
sy pi fe ab]? Vv 
2 2; 2 2°) 
t bTF ] 3 V/ 
ames Lda 
But 
j a b l 
= — Mt DD » abd, 
D,* +” M, — D,D,D, —| DabM 
| MM, = DD, + D,D,+*t” m,, 
lp ow, D,D.D, — D2D,— D,2 +” M,. 
by substituting these, the equation becomes 
: ata b l ‘ dia b } 
D2D.—D 1 3)" ab] M, {1 4 ab ah | 


] j 3M? a brs 3M 
| i v, | aos | 8H | 0, 


This may easily be transformed into 





D2D, D1 D? + 8D,D, ~— b p21] 
DD:| : = ab 


») 


p,| De 2D. D, 


aes "ls D? +3D,D,—2 YP!) = 0. 


s 2 MM, 


~) 
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Which reduees to 


8a+bDD, , 3D7—1 5 4 »§ @ ery 
[3 22 M2 M, 2 Al . i | 
a br3 ; 3 D?D-? 

D9 ft 7 | *@, 


and thence to 


] 
or 


which is perceived to be identical. 


D?—1 3 
D?+1 2D7+1 


bo Ce 
~~ 
~~ 
Cansesseeal 
~ 
_’ 
~~ 
PE 


i a 


W hen 
l ‘ 
kh, Fva i 
| ( 
is divided by 1] v the remainder is equivalent to what 7 becomes when wz 
1 & 
is put equal to 1. sut 
] (a rite b) 


] ( \ ] rr HT? 
consequently this remainder is 


(1 a) (1 b) 


/1(1) 
the ambiguous sign being so taken as to render the quantity positive. In like 
manner it is shown that the remainder of 7 divided by 1 +. is 
J 


1 (1+a)(1 b) 


Then 


1d (@ y Ct) 1; (0 a) (1 b) 2 
dd = : 1 : . af ay 


where the upper or lower sign is to be taken according as: //(1) is positive or 


negative. And 


d (a ; Ct) 1; (1 a) (1 b) I’ 
dy 2 ] L 





@ 


ee 
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where the upper or lower sign is to be taken according as //(— 1) is positive 
or negative. The expressions for the / and /’, correct to quantities of the 
order of the fourth power of the eccentricities inclusive, are 


we dt), 


aK. + Dot Dr DF HAE. 


Lei. D.|h, [arn D, D|E IT) E., 


[oA D.—2Dy |B + [RGR — De] B+ Dube 


a ee 
A ah 


i dD, . r dD, ; 


| A 7 D, 20, | Ey |“ oe D,| E, DE, « 


D, + D, D,| H(—\(E 





By substituting the value 


and putting 








; 
| 
' 
| 
; 
: 
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where we have, as has been proved above, the relation J, N,, we get 
So . £7) i 44 ay] ly) r ’ ’ / 
. ¥ N NV, cos AN »>COS DN 
dd a b a 1) P - ° é 
: | ) COS C’ 
dia L C2) Ly (1 a)il b) . : ; “ ? 
_* . N, NJ cos ¢ N,’ cos 2¢ 
Ad at b il b ‘ Y, 
. ] 9 5 COS ¢' 
Integrating, we have 
’ | a 7 rag sie a pets , 
wo te. ( are tan tan? Nid N, sin ¢ 1 NV,sin 2¢ 
’ rN | by Z , 7 
' y [ l a dn] Lo iV . / ' V — - 
o Lct—=. are tan tan ’ N,'¢ NV sin ¢ LV,’ sin 2¢ 
: Niyb*"9 pa 








st se 


Sa°8r 


The quadrant in which the are correspondent to the tangent is to be taken 


is found by dividing the number of the quadrant of ¢/ by 2, if it is even; 


or 


by augmenting the number of the quadrant of ¢) by unity, if it is odd, and 


then dividing by 2. 


By taking the sine, we have, 4 being any arbitrary angle, 


1 | rsin(@ 1% 4+ 3) 
. a rs : ' >» yy 
» | asin 5 cos| Vc + ¢ + 3+ sing VV, sin 2 
o.. - re , , - 
1 1—beos5sin[ Mi +e | 3+ sing + £2, sin 24 
| ] rPSIn ” Ct 7) 
/, 
4 —— a a — ae 
1 1 sin 5 cos Nv \3 NY sin 'N,’ sin 2. 


/ 


 . v, ‘ — as : 
bcos} sin[ NV,'¢ + ¢ + 3+ sin ¢ +- ,X,'sin 2 








VA eee 
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or, as they may be written, 
; | “SID (@ yO j3) 
iy 1—bs4 1—alsin[(( V4) e+ 35 Msing 4 4N,sin 2% 
ify 1—b+) 1—a)]sin[(V,—))¢ e434 V sing + 4N,sin % | 
yy] rsin(@ LC ¢ -+- 3) 
at l b cg a}sin | N x4 ( 3 Vi sing : V, sin 2¢ 
Ty lt+b+) 1 a] sin | V fh --e'+-3-+-D,' sing +3 1,’sin2 
The expression for the auxiliary angle ¢’ in terms of the time, which has 


tlready been obtained, we will denote is follows 


ru f(t c.) A, sin #, (¢ ci) A’, sin 24, (¢ 


Substituting this for ¢ in the preceding formule, and putting in succession 


4 Sha J ) LL 
We vet 
sin sin 
| ] r D) tyod L | | i /? -)O,(0 
COS ° COS ” 
JP? sin 4, (¢-- P. sin 24, (7 
sin 
1} 1 ly y J a ar W(t 
; cos ’ 
: f \ ein 8H if 
(/ sin #.() ( in 2 
sin . sin 
] l r w std | ly ] | a | [(P at ( 
COS a COS , 
IP’ sin f(t «,) Py sin 24 (¢ 


a} (UP | Af a\ 


COS 


(), sin“ (¢ 


; ) j 
(/ sin vl 
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Here we have put 
N, 
P N, + (N+ 4)4,, 
ae b| .V, + VA, + 2(N, + 3) AL 
(, N, (N,, 1) K,, 


Y, = 4[N,4+ VA, + 2N,—)A, 
Pi=N/4+(%4+)4,, 

Py = }4[N, + WK, + 2(N,4+ 4) All, 
Y/= N+ (N,— pb), 


cl | ( 


Q/ =i Ny + N/K, + WN, — 4) K. 


It is evident from the equivalent of sin » cos 7 derived from these equa- 
tions that e eore C 180°, according as 


TT (hb) D), Db DIP Dp ieee j l /- 


is positive or negative. Hence the latter of the two equations may be written 


j sin . : sin D i ; : 

| | ie i) s[y l b ] l a | oan [ Fr. 3) aii¢ Ga ( 

Pi‘ sind,(¢-«,)- LP, sin 24, (t- ¢,) a 
i _7 sin . ; 

ify l+b+) 1 ®] og Lie S)A(t- e)- 

(),'sin 4, (¢— ¢) 4+ @,’sin 26, (t+ ¢,) 4-...], 
; 

. 


where the upper or lower of the newly introduced ambiguous signs is taken 


according as //(b) is positive or negative. 


Let us put 
(Po +4 
x =(F,.—4)4( +4) +6, 


™ 


| 
» 1—biyl a], 


+ 
[ 
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Then 
sin sin 
] i ” ( | 72) () 
3 cos - t/ 1 lel! eos 
4 sin 
: (1 oe 1 : 
Fe l 4¢l | cos 
sin 
1 > ] ] ) . 
LAP, (AY, L(,) ues 7) 
¥ »? ) sin 
“9 | 4 iV, (4 1 L/ , 2 /) 
a COS 
4 ] , ] , sin ) ) 
, sin ©. ) 
( / 4 (/ ) eon (97 =/ 
sin si 
os 1 7D) (1 LP?) } i) — 
COs ‘ , - . cos 
(] 1¢)”) | P Sin 
COS 
LAP, L¢) 1yy’)| 5M 9 
2 . vltiihas COS 4 
LAG P 1P,)/7" (2 / 
~ ilies a ‘a COS 
1D Lp sl (2 Day’) 
(42? =" =/ eos /. a 
] ] SIM (a. +). 
{ (), 5) | etn \°R =/ 
: ° ° sin - sin . ° ° . 
4 It is evident that « o and « wo can be expressed in series of the 
a COS COS 
same form. 
In applying to Jupiter and Saturn these equations, it is found that by 
varying the value of A, 
2s -d AD . AD 
a re D, 0.0101629, A “% 0.0009178,  Av' / 2 0.0050568. 
: dh dh dh 
Also 
log | ; (1 a)(1 b)| 9.9574334n, log) (1 a) (1 b) Q 4074864. 
a Ll, 0.1672972, ; 0.1655301, 
L, 0.0028107, L, 0.0012760, 





0.0000071. | i 0.0000250. 




































a a ae 


W hence 


Also 








rHE 


0.16676382, 
0.0021649. 


0.000002 1, 


0.6837293, 


9 5750158, 
O.0OL01L623. 


99" 55981. 


The value of ¢ is found to be 


Hence the expressions for the two arguments are 


0.3759635 °" ; 
COS 

0.0018085 2" 
COS 
0.0000072 72 
COS 
0.72938089 7 
COS 
0.0013889 72 
COS 

0.0000058 EY: 
COS 

e — 340° 


308° 13’ 15”. 


31 ; ae 
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LD. 


() 
‘ 


(/, 
v2 


log 


log 


(ft; 


+). 


8’ 50 


7) 


“7 ) 


7) 


) 


27) 


26. 





0.16676382., 
0.0009746, 


0.Q0000074. 


9 8634412 


9. 7217366 


L) a. 3° 5OLS6. 


2 


1.024ys24 SI, 


COS 


0.0001196 2 


COS 


0.0000016 22 


COS 


0.5255160 


COS 


0.0008842 S12 
COS 


nh 


0.0000052 7" (: 


COS 


22” 5598172, 


3”. 501562. 


(27 


7) 


) 


27). 


47), 


oes 
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The following expressions for ¢ and e were obtained : 
t 


; 8.6282138| , (1 — [6.5410419| cos ¢), 
] 7 
8.8231642| , (1 + [6.9312571 | cos «), 
yi-+a 
r 8.6282135 | 16.2400] cos (7¥—7) 3.7900 cos 2 7 7) 
—— 
, S.8231648 | |] 6.65023 | cos (7 7) 1.1982 | cos 2 (7 7 : 
as 


By means of these we can pass to the expressions for the following functions 


Sill F 0.01596822 Sin - 0.04354278 Sin 
COS Coe °* COS 
~ (0¢ , Sin °. “y’\ , a SID 2. f 
0.00005696 7° (27 / 0.00000886 "(27 / 
sin s1n 
0.00000031 (33 By’) O.Q0QQ000009 (97 2. 
COs 7 VA COS ‘i “f)> 
sin . - sin yy SID 
» 0.04852990 . 0.0349640 , 
eos: , COS A i 4 , COS A 
- sin - sin 
0.00008205 ~" (2- ~7') 0.0000513 2. 
’ cos A A I 4 COS az VA 
— ’ © IM o 
0.00000033 (87 — 27’) — 0.00000031 * , (By 27). 


cos COs 


It will be observed that these expressions are as convergent as could be 
wished. The form of these integrals being discovered, another and mor 
direct method of arriving at them is suggested. The coefficients being 
assumed as indeterminate as well as the rates of movement of the two 
arguments together with the constants which complete the values of the 
latter, the expressions could be substituted in the differential equations, and 
thus would arise twelve equations of condition, which along with the value 
of the four variables at the origin of time would determine the sixteen 
unknowns involved. But on trial it seems this way of proceeding would 
necessitate as long computations as the method we have followed. 

In conclusion, it may be observed that, if terms arising from the squares 
and higher powers of the masses were taken into consideration, the form of 
this investigation would not thereby be changed; the only effect produced 
would be that the values of the various constants involved would receive 


slight modifications. 








a 
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About the vertices of an equilateral triangle three spheres are drawn with 
radii equal to the side of the triangle. Find the volume common to them all. 
W. WM. Thornton. 


SOLUTION. 


Let ABC be the triangle of the centres, V its ortho-centre, and V one of 
the apices of the solid. Draw LO to meet the surface in Y The planes 
VOA, NOY, AVY, and the sphere-surface AVY bound one-twelfth of the 
solid. The area of the sphere-surface A NY is the ditterence between the area 
of the zone-surface 1¢(’V, whose angle A CN is are tan 24 2, and that of the 
right spherical triangle YC, whose angles are are tan 2) 2 and are cot 2) 2. 
The excess of this triangle is 
/; = are tan 2) 2 are cot; 2 a i 4 
Its area is 4 
a ER? 
/? being the radius of the sphere. The area of the zone-piece is 
& » P? are tan 24 2. 
Hence the curved surface of the portion of the solid under consideration is 
Z—dJ I? (4 =z — are cot 1 2 y are tan 2) 2) : 


/?* (are tan 1 2 S are tan 2) 2). 
The volume of the spherical pyramid which has this surface for its base is 
\ /2* (are tan 1 2 Sare tan 2) 2). 
The volume of the cone-segment whose base is VA and apex B is 
: I? f are tan 24 2 v7 » 2). 
The difference between these volumes is one-twelfth the volume of the solid 


common to three spheres. Hence 
V = f(s) 2 + 4 are tan) 2 — 3 arc tan 2) 2). 


W. HL. Echols. | 


Ga ic 


4 a 
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a h, a. < r Ab 

a, b, c eo h, <a ag — 

au, h, Cs sas h, a, h a, ¢ a j 

dj =| Gy h,, ( ee h., , and / U. h a, ¢ 1, I 
se & & h,, a, b a, ¢ - ah 

a, > ad, € rh 


then will dD a” “2. TZ. M. Bluksle 
SOLUTION. 


This exercise is given on p. 77 of Muir’s Theory of Determinants. It may 
be got from the result of section 55 of that work; or it may be derived by 
a process analogous to that of the section referred to, as follows: Multiply 


each column after the first by 4,; add to each element of the second column, 


thus multiplied, }, times the corresponding element of the first column; to 
each element of the new third column c, times the corresponding element 
of the first column; ... to each element of the new ath column 4A times the 


corresponding element of the first column ; and we have 


i, ) () () 
“/, ah, a, ae udoe = ud hj ad /, 
a,” Jaz! a, ab, a,b a,€ At i ah a.h 

a, ad, h, a,b L,C, d¢ = if /; 4.4 

uf /, a, ( ‘ if} 

(l, /, a, Ce a “ h. 
(ty aD; 

a, h, , Cn —o a, h,, 

“) J PD. 


W B. Li hards. 


Also solved by L. G. Weld and the proposer. 
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Show that the attraction of a finite mass on one of its points is finite. 


° A. Hall. 
SOLUTION. 


Take the origin at the point attracted, and let /drdydz be the element 


of mass, |’ the potential for a particle of the body, and 7 the distance of the 


particle from the origin. We have 


bd bd hdr dy dz 


foe: | i ieee 


Put 
/ 7 COS % , y/ ? Sin 7 COSA, rsin sins; 
then 
dr dy dz ia sin uf du di d) . 
and 
. . . 
| Ker sin uu du d) di ‘ 
The limits of integration are 0. to z: 4=0, tos 9-7: and 7 0. 


to the limits of the body. For the component of the foree in the axis of « 


we have 
° |” e Pp hr dxrdy dz rf pf 2 . ’ 
x : j ‘OS Sin 2% du didp. 
- “Seo fy freo 
with similar values for }” and Z These components are finite, and therefore 


the le sultant is tinite. See (vauss. A llgeine ine Lehrsdtze “ ere. 
A. Ilall. 


EXERCISES. 
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A body at distance 7 from the sun is moving with velocity 7. Prove that 
the major axis of the orbit described is parallel to the direction of motion if, 
and only if, the velocity is “circular velocity for the distance 7.” 


Ellery W. Davis. 
307 


If a horizontal beam of length 2a is supported at each end, and has a 


load in the form of an isosceles triangle, base 2a, height 4, a unit’s thickness 


throughout, and heaviness unity ; show that the deflection of the beam due to 
2a'h 


Lk] [7 i, Luylor. 


this triangular load is 





ee 
“ae 
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